iL9,3L0 


NASA-CR-169360 

19820025547 


A Reproduced Copy 

OF 


Reproduced for NASA 
by the 

NASA Scientific and Technical Information Facility 


LIBRARY GSPY 


1384 


FFNo 672 Aug 65 


UNGLEY RESEARCH CENTER 
LIBRARY. NASA 
HAMPTON, VIRGINIA 



\ 






HOWARD UNIVERSITY 
SCHOOL OF ENGINEERING 
DEPARTMENT OF MECHANICAL ENGINEERING 
WASHINGTON, D,C. 20059 


FINAL REPORT 

NASA GRANT; NSG— 1414, Suppl, 4 

THE DYNAMICS AND CONTROL OF LARGE 
FLEXIBLE SPACE STRUCTURES-V 



by 

Peter M. Bainum ^ 
Professor of Aerospace Engineering 
Principal Investigator 


and 


A*S«S*R« Reddy 
R* Krishna 
Cheick M. Diarra 
V*K. Kuinar 

Graduate Research Assistants 


August 1982 


I 













ABSTRACT 

A general survey of the progress made in the current and past 
grant years In the areas of mathematical modelling of the system 
dynamics, structural analysis, development of control algorithms, 
and simulation of environmental disturbances is presented. The use 
of graph theory techniques is employed to examine the effects of in- 
herent damping associated with LSST systems on the number and locations 
of the required control actuators. The presence of damping allous a 
greater flexibility to the selection of actuator locations under which 
the system is controllable, while the rank characteristics of the system 
influence both the number and locations of the required actuators . 
A mathematical model of the forces and moments Induced on a flexible 
orbiting beam due to solar radiation pressure Is developed and typical 
steady-state open-loop responses obtained for the case when rotations 
and vibrations are limited to occur within the orbit plane. A preliminary 
controls analysis based on a truncated (13 mode) finite element model of 
the 122m. Hoop/Column antenna indicates that a minimum of six appro- 
priately placed actuators is required for controllability. An algorithm 
to evaluate the coefficients which describe coupling between the rigid 
rotational and flexible nodes and also intra-raodal coupling has been 
developed and numerical evaluation based on the finite element model of 
Hoop/Column system is currently in progress. 
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I. INTRODUCTION 


The present grant represents a further extension of the effort 
initiated in previous grant years (May 1977 - May 1981) and reported 
in Refs. 1-7*. Techniques for controlling both the shape and orienta- 
tion of very large inherently flexible proposed future spacecraft 
systems are being studied. Possible applications of such large struc- 
tures in orbit include: large scale communications; earth observation 
and resource sensing systems; orbitally based electronic mail trans- 
mission; and as orbital platforms for the collection of solar energy 
and transmission (via microwave) to earth based receivers. 

This report is subdivided into seven chapters. Chapter II is 
based on an invited general survey paper presented at the recent 10th 
IMACS World Congress on System Simulation and Scientific Computation, 
-August 1982, and presents a general survey of the progress to date in 
four general areas r (1) mathematical modelling of the system dynamics; 
(2) structural analysis; (3) development of control algorithms, and (4) 
review of previous work in the simulation of environmental disturbances 
(mainly due to solar radiation pressure) . 

In Chapter III, the use of graph theoretic techniques, previously 
introduced^ to simplify the eigenvalue calculation for LSST systems by 
reducing the system matrix to a collection of lower order sub-matrices. 


*For references cited in this report, please see list of references 
at the end of each chapter. 
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is extended here to address the controllability of . inherently damped 
large flexible space systems. A second paper to be presented at the 
33rd International Astronautical Congress, Sept. 1982> forms the basis 
of this chapter. 

At the operational altitudes of the future missions involving large 
space structures, the principal environmental disturbance is that due 
to solar radiation pressure. The effect of solar radiation (pressure) 
disturbance on a flexible orbiting free-free beam is addressed in Chapter 
IV, and to the authors* knowledge represents a first attempt to include 
such disturbances in the system dynamics a flexible structure in 
orbit. (A paper based on Chapter IV has just been accepted for presenta- 
tion at the 1983 AIAA Aerospace Sciences Meeting, January 1983.) 

8 

Our proposal for the 1981-82 grant year originally emphasized work 
to be performed in three areas: (1) further analysis of environmental 

effects; (2) graph theory approach to the controllability, observability 
and eigenvalues of large scale systems; and (3) consideration of sensor 
and actuator dynamics. Shortly after submission of this document, we 
were advised by NASA-LRC that it was desired to redirect our effort so 
as to provide direct support to the synthesis of control laws for the 
LSST Hoop/Column Maypole Antenna system whose feasibility is currently 
being studied by the Harris Corporation, Melbourne, FloriV'. As a result 
of this the third task listed in our proposal was not considered, but, 
instead. Chapters V and VI represent our preliminary efforts in support 
of the Hoop/Column controls analysis. 


r; ! 
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in connection with this the OhAas conputet elsotlth.’, developed 

by KASA LEC Is being used entenslyely In this effort. A •■VersetUe Hoop/ 
C010« Antenne Sttnctersl Hynenlcs Hodel“" based on . HASIRAN finite 
Clement software package, prepared by the Harris Corporation, has been 
transmitted to us by HASA-IRC. togehter with a magnetic computer tape 
containing the eigenvectors for the first 34 modes of a single layer 
surface model of the 122m. s»del of the Uoop/Column - the latter received 
shortly before the end of the 1981-82 grant year. The effort described 
In Chapters V and VI Is being continued during the 1982-83 grant year. 

Chapter VII describes the main general conclusions together with 
recommendations for further work. 
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II. ON THE MODELLING AND SIMULATION OF THE DYNAMICS AND CONTROL OF LARGE 
FLEXIBLE ORBITING SYSTEMS 


This paper attempts to review the steps Involved In the development of mathematical 
models that can be used to simulate the In^rblt dynamics of large flexible systems. 
The use of graph thuicetlc technl<{uas can often be used to reduce the computational 
effort Involved for calculating the eigenvalues of large ordered systems. Computer 
generated graphical techniques may provide additional Insight Into the understanding 
of elastic ssdal shape functions of complex systems. Finally the nuser leal techni- 
ques commonly used to develop shape and attltxide control lavs vUl be briafly reviewed. 
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I. INTRODUCTION 

Large, flexible orbiting systems have been pro- 
posed for poselble use in communications, elec- 
tronic orbital based mail systems, and in solar 
energy collection. The size and low weight to 
area ratio of such systems indicate that system 
fle.xlbllity is now the main consideration la 
the dynamics and control problem as compared to 
the inherently rigid nature of earlier space- 
craft systems. For such large flexible aystems 
both orientation and surface shape control will 
often be required. 

Fig. 1 illustrates a conceptual plan of develop- 
ment of a system software capability far Use in 
the analysis of the dynamics and control of 
large space structures technology (L5ST> systems. 
This concept can be subdivided into four dif- 
ferent stages: (1) system dynamics; (2) struc- 

tural dynamics; (3) sppl lest ion of control algo- 
rithms; and (4) the slsrilatlon of the environ- 
mental disturbances. The most fundamental 
ponenc Is that of the modelling af the system 
dynamics of such systems in orbit. 

II. MATHEMATICAL MODELLING OF SYSTEM DYNAMICS 

Previously many authors analyzed spacecraft sys- 
tems coosi^clng of a primary rigid body and 
elastic appetlages which represented solar pa- 
nels, antenna booms, last rumencac ion platforms, 
etc. The hybrid coordinate Dodelllnc; method as 
Incroductcd by Mairovitch^ and Nelson^ and fur- 
ther developed by Lilcins^ has been widely used 
In the study of such systems. The hybrid coor- 
dinates are comprised of the attitude (Euler) 
angles or quasi-coordinates of the primary body, 
as discrete coordinates, together with the dis- 
placement of the elastic appendages relative ca 
the primary, characterized as distributed or 
iDodel coordinates. 


^Research supported by NASA Grant; NSC-L414 
*^Professor of Aerospace Engineering 
**Graduate Research Assistants 


The presence of rocers and movable damping 
mechanisms on the main pact can be readily in- 
corporated and the error involved in using a 
finite truncated series representation decreases 
with an increasing number of mode shapes. The 
remaining high frequency errors should be atten- 
uated by suitable design of Che closed-loop 
attitude control system. 

Vhen the siza end vrels^t to area ratio of sxich 
proposed future LSST systems indicates that the 
entire system oust be esasidered to be flexible, 
the hybrid coordinate formulae ion laay not always 
be readily adapted by simply assuiJlng the mass 
and Inertias of the (previously) rigid central 
part tend to zero in the Unit. Santinl^ has 
developed a sachemtlcal formulation for pre- 
dicting the motion af a general orbiting flexi- 
ble body using a coatlmum appi oach. Elastic 
deforeations are consldtrad small as compared 
with characteristic body dimensions. Equations 
are developed for both the rigid and elastic 
(generic) modes. This development is based os 
an a priori knowledge <ef the frequencies and 
modal shape functions <of all oades to be in- 
cluded in the truncated system model. Kumar** 
and Bainum** have modified the development of 
Ref. 3 based on vector calculus so that elastic 
modal shape functions ^expressed in arbitrary 
systems of coordinates say be accommodated. 

III. STRUCTURAL ANALYSTS 

For simply isotropic structures, such ar homo- 
geneous beams and circular places, closed form 
expressions are available for cha elastic modal 
frequencies and shape i^snetions. For more com- 
plex and/or noalsotroplc sy items numerical 
methods must be employe! to obtain this informa- 
tion. Commonly used roacines include versions 
of STRUDL and NASTRAN Ixhe latter more complex 
algorlcha cay also be ixceful In Che simulation 
of chertDclasClc envlrcrsental effects). The 
use of computer generalsd graphics may also 
prove useful In understanding the elastic shape 
fmctloas of complex sfscems when excited at 
different codal frcqucstlcs. 
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Fig. 2 UXotcrjicis chs sod&l ^hap# funccloas of 
1 fr«e-f^-i« circular alaaclc placa, showing ths 
prssssca of ths nodal linas and rMrldlans as a 
fuacclon of cha frequency peraaecer, X. 

For syscens requiring a large number of elastic 
nodes to accurately represent the systea dyna- 
alcs, graoh theoretic techniques may provide an 
alternative to the numerical problems Involved 
In calcxilatlng the eigenvalues (modal fre- 
quencies).’ With this approach the system 
(stiffness) matrlJC can often be reduced to a 
system of lower ordered subcacrlcea so that, 
under certain conditions* the eigenvalues of the 
original matrix are given by the union of the 
eigenvalues of the subnatrlces. Including their 
culclpliclty. For large order systems this 
approach can substantially reduce the numerical 
effort Involved with an ieprovement In accuracy. 
As an example* a free-free homogeneous 
square placa was considered in Ref. 5 where 
Che original stiffness matrix contained a dimen- 
sionality of 16. Although this matrix was 
sparse (*ess than 20Z of the elements 
were non-rero), many of the non-zero elements 
were off-dlagonal. Fig. 3 represents the di- 
graph of cha 16x16 stiffness matrix and Indi- 
csces chat, under appropriate conditions, the 
eigenvalues of the original matrix cay be ob- 
tained by calculating the eigenvalues of seven 
sub-iatrlces; the largest dlnenalonellcy of any 
of the sub-eacrlces here is 4x4. Complete de- 
tails describing this example and an algorlthr* 
that may be used to determine the sub-matrlcos* 
if Che original matrix Is reducible, are pro- 
vided In Ref, 5. The non-zero elements or cho 
■stiffness matrix are given In Table I, The 
eigenvalues are evaluated from the original 
16x16 stiffness matrix as well as the seven 
reduced order sub— matrices and are compared In 
Table II. 


TABLE I 

54jn-zero elements of the iC-nacrlx: 


iC(3,3)- 34.919 
K(12,4)« -19.9054 
K(14,6)« 15.4658 
K(15,7)- -11.6130 
K(16,8)« 96.9152 
K(4,10)- 17.5471 
K(9,1D- -5.31786 
K(12,12)« 326.604 
K(6,14)- 15.4658 
K(7,15)- -11.6130 
K(3,16)- 96.9152 
K(ll*3)- -5.31786 
K(6,6>- 13.3938 
K{7,7)- 34.9190 
K(3,8)- 336.913 
K(9*9>- 34.9190 
X(12,10)- -11.6130 
K(4,12)- - 19.9054 

K(13,13)- 265.337 
K(14,14)- 386.913 
X(15,15)- 325.604 
K(16,16)- 1732. S5 


K(9.3)- 8.03522 
mo. 4)- 17.5471 
K(16,6)- 17.8584 
K(6*8)- 15.4658 
K(3,9)- 8.03522 
K(10,10)« 34.9190 
mi.lD- 76.8770 
K(7,13)- 17.5471 
K(8.14)- 56.1775 
m3.15). -19.9054 
m4,16)- 96.9152 
K(4.4)- 265.337 
K(8*6)- 15.4658 
m3, 7)- 17.5471 
K(14,8)- 56.1775 
K(ll*9)- -5.31786 
K(3. ID- -5.31786 
K(10,12)- -11.6130 
ms. 13)- -19.9054 
K(16.14). 96.9152 
K(6.16)- 17.3584 


TABLE II 

Eigenvalues of the K satrlx: 


Submatrlcfa Original catrlx 


S, 

J.O 

0.0 


0.0 

0.0 

*3 

26.88373000 

26.38378620 


41.36168802 

41.36169439 


78.46953198 

78.46953100 


33.24206125 

33.24206653 


260.2368175 

260.2367428 


333.3811212 

333.3811150 

*5 

0.0 

0.0 

*6 

12.21269086 

•2.2126»085 


330.7355000 

jJO. 7355804 


429.6358695 

429.6359588 


1747.485740 

1747.485839 


33.24206125 

33.24206653 • 


260.2368175 

260.2367428 


333.38U212 

333.3811150 


The eigenvalues of the system are the eigen- 
values of ir^K and as .M-ml, elg<i»*'/alues of 
.'T^IC-d/n) times eigenvalues of K. 

Before surface and orientation control systems 
can be designed, It Is necessary to understand 
the dynamics and stability of the 'mcontrolled 
system. For Virge order systems an analytical 
approach to the stability problem is not fea- 
sible and numerical techniques must be employed 
to develop the system characteristic equation 
and the loci of Its root* for different sets 
of system parameters.* As the number o' modes 
retained In the truncated system esdel In- 
creases, expansion of ch« characteristic deter- 
nlnental equation becomes algebraically pro- 
hibitive. As an alternative an algorithm due 
to Leverrler • • can be used to numerically 
determine the coefflcler.es tn the charac- 
teristic equation. In erde;.- to Implement chls 
algorithm the linearized equations onist be 
written In standard state variable format, 

IV, CONTROL ALGORITHMS 


AC cais point the trodelllng of the control 
actuators can be added to the pr^^vloualy de- 
veloped open-loop system codels.**^^ In gen- 
eral, an actuator placed at an arbitrary loca- 
tion on a large space structure will affect 
both the rigid and flexible ojdes. The loca- 
tion of auch an accuacor has daiflnice lapll- 
catlons on ch- system controllability. Fcr 
large ordev systems the reachabUlty* matrix 
and term i-ank concepts, also developed from 
graph theoretic techniques, may be used to 

can be loiouta- 

tionally more ef..ectlve chan numerical rank 
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ccists of the system conctoXlsbiXlty matrix. 

Tim. free-free square piece la esaln ccaaldered 
es en Uluetretlve e x a mp le In Eef • 5. 

Tlxree cecfaalques ere cornnnly used to develop 
control levs once the syscsM cancrollebllicy 
has been esteblished. These laclude: (e) de- 

coupling techaiquee; (b) pole pUccaeat (clus- 
tering}; end (c) ea eppliceclon of the lineer 
regulator pcoblca from optimal control theory. 
These three tecfcalques approach the controls 
problem froe different polnte of vlev end each 
vlll be briefly discussed. 

The decoupling techalqut cen be applied In two 
distinct sub-cetes: (1) where the linear state 

equations la the original coordinates ere de- 
coupled by using state variable feedback tech- 
aiquee; and (2) where the open loop linear 
equatlone ere first treasforaed Into e de- 
coupled set la modal coordinates and then con- 
trol laws are developed Independently for each 
m>de. It then becomes necessary to transform 
the control lavs as expressed In modal coordi- 
nates to the actual control In the original 
coordinates. 

As an example of sub«-case (!)• we assume that 
the llaeariced equations can be expressed as: 

Z - DZ + E2 + BU (1) 

where Z • ****r+n^ 

describes the rigid body position dlspXacercnts 
(1,2. ..r) plus the a elastic (position ) 
coordinates retained, la any truncated codol. 
After selecting U-K^-HCpZ, wo can rewrite the 
controlled ODtion eqiiatloas as: 

Z - ^ (E+Brp)2 U> 

where Z and t, are evaluated such that (IHBTCf) 
end (E+5 jU) erS dtaEonalUcd and thus yield the 
required canplng and frequency of the con- 
trolled osdes. The total number of rlgid+ 
elastic modes, (r+n), must be equal to the 
number of actuators here to avoid the necessity 
of using pseudoinverse matrices. 

For the second tub-case, as an example, let us 
consider a different fora of Eq. (1) In terns 
cf the mass (M) and stiffness 00 matrices, 
where D-0. 

M r rz • F - BO (3) 

Vlch the following type of transformation: 

2*0q, Eq. (3) say be recast la terms of the 
modal coordinates, q, and the transformation 
matrix (Involving the eigenvectors). 0, as: 

such chat and are diagonallxad 

C'MjJ q f [‘Cj^Jq - - r (5) 


It is thiKi possible to design the con*:tol laws, 
F*, In the modal space so that ladependeat con- 
trol of each of the modes can be achieved. A 
transforcatlon is then required to obtain the 
control laws la the original coordinates, F, 
and, then, for a given location of ecti^cors, 
the actual control from: 0-B*^F-B*^(^^)"‘^F*. 

In the pole clustering method the overall tran- 
sient requirements if the system are considered 
Instead of coacentraclng on the behavior of the 
individual coordlaatiis. The linear used aystea 
equations, Eq. (1), ^:att be recast In the state 
space format as: 

X . AX+SU W 


where 


m • • • • 


The control, IK-I3C, is then selected by uslnj 
a digital computer algorichu such as OBACLS* 
such that (A-BE) has the identical negative 
real part In each of Its eigenvalues. Al- 
ttnugh the number of actuators can be less than 
the nu^er of modes, e limitation of this pax- 
tlcular algorltha is that the gains are se- 
lected such that all of the cloied-loop poles 
lie on a lino parallel to the Imaginary axis. 
The alcoritha is useful, however, when it is 
important that each mode in the system satisfy 
some damping characteristics. 


The linear regulator theory allows the analyst 
to sat, a priori, distinct penalty wel-thr'-ng 
functions on the control effort as well as the 
state variables. The control law, C* — CC, is 
salecced such that the following performnee 
Index is mlniaizcd 

j - ;(?Qx u^atndt V) 

o 

where Q and R are positive definite penalty ma- 
trices. The steady state solution of ths ma- 
trix Rlccatl equation of diaenslon equal to r.he 
state has to be solved In order to obtain the 
gain matrlcss, K (which represent the positive 
definite soLution to the algebraic matrix 
Riccatl equations). 

A computer algorltha within the CfLACLS^* soft- 
ware package can be used to obtain the gain 
Bat rices, K, for different combinations of the 
Q and R penalty matrices. This algorltha uti- 
lises the Mewton-Raphson method of solving the 
Rlecatl equation. 

Both the linear regulator problen and the pole 
clustering methods can result in sqm of the 
closed loop frequencioa being orders of sagai- 
tude greater th^ those of the uncontrolled 
syscea. T^ee« hlgh«r £r«qucocl*« al»o 

correspond to the frequencies of higher modes 
not Included in the previously truncated sys- 
tem 3 »del. In order to cotrpltctly consider 
such effects the order of the original system 
sodel would have to be Increased in order to 
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tbtt of coatrol tplllovor. On cho 

othor hsnd» choto eocbodi hav« eho advantaso 
thac choy car ba applied to aUuaclsnt vharo 
Cha nus2>ar of actuatora It lota chan cha nua- 
bar of oodaa In cha sachceacical aodal» in 
concraac to cha uaual applicaclona of cha da* 
coupllns aachoda. £xa8:piaa of cha appllcaclon 
of cha aarloua control alcoclchssa ara $lvan in 
Kaf. 10. Aj an axanpla, a typical application 
of daccupllne, ualng acaca varlabla faadback 
for cha orlcacacioQ and ahapa concrol of a 
a frcffraa aquara place nooinally follovlnf 
cha local varclcal with lea larger aurfaca 
noraal co cha orblc norsal la conaldarad. Tba 
eaodal eoncalna thraa rigid rocac tonal eodea and 
Cha first chraa cranrrarsa flaxlbla sodas, wlch 
alx accuacara asausad ca ba located at shown In 
Fig. 4. The decoupling galna ara salaccad in 
order to peoduca of cricical daaplng In 
each of cha rigid oodaa and cha fundacancal 
alaacic node and lOi of critical danplng In 
Cha aacono and third flaxlbla oodaa. Tha 
cone roll ad scat a retponta la given In Fig. S 
and cha corraspondlng cioa hlscory of cha ra* 
quirtd concrol forcaa la llluscracad in Fig. 6. 

V. SIMULAIIOS OF EhTIROSXErTAL OISTURSAXCFS 

Tha principal dlscurbaoca forcaa and corquas 
acting on a larga flaxlbla ayscem In orbit ara 
Cha gyroscopic and gravity grad lane corquas 
assoclaced with cha orbital aotlon. cha concrol 
Corquts, and chosa corquas due co cha tnvicon- 
teic. In cha forsulaclon of Raft. 3 and 4 cha 
gyroscopic and gravity grad lane corquea ara 
includad In cha xodal of Cha ayaceia dynaiilcs. 

If othar formulae Ions an aaployed. such as 
gtnaral flnlca alasenc atchods. which do not 
account for cha orblc&l dynaslcs. cha effaces 
of cha gyroscopic and gravity grad lent corqiias 
should ba carefully considered before delating 
Chen froa cha dynaalc nodal. Tha treacs^c of 
concrol aooalllng and algorlchas was axaalned 
in cha last sect Ion of chls paper and no fur> 
Char alabo ration will ba provided hart. 

EnvlsonBancal disturbances can ba accrlbucad 
aaisly co cha affects of solar radiation pras* 
sure, except In very low earth orbit where tha 
aarocyna&ic drag forces prcdoalnata. MoDancs 
due to aolar radiation prasturc ara Induced If 
tha center of solar radiation preasure is not 
co<-Iocac«d wlch cha syscea center of saaa. 

The location of the center of pressure la de- 
pendent on tha surface charactarlitlcs as well 
as cha gaocecrlcal shape of the structure. In 
addition, due to solar heating, ther&al gradi- 
ents can ba induced In Cha structure which say 
result in appreciable chersnl strains. As a 
result, the struccura will undergo deforaat lor.s 
which will further concrlbuca to cha forces 
and torques caused by cha solar prasaura. 

Several Invest Igacors have considered the af- 
fect of solar radiation pressure on the dyna- 
aics of spacecraft. The aajorlty of the space- 
craft toddled conaisetd of a aaallar rigid 


ceacral secaUlca co which flat plate append- 
ages. also craacad as rigid, were aasusod co be 
attached. A few aushors shoved how cho eolar 
praseure conents geaerac^ could be used for 
aacelllte attitude control by concroLllng Che 
orientation of places and/or vanas which could 
rotate at the ande of the appendagei^* * ’ An 
axtanaloa of chaea codeia ca Includa large In- 
herencly flexible orbiclng syscecs la naedad 
before the nacura of cha snvironeeacal dls- 
curbancaa on propoaad LSST ayscaas can ba coa- 
placaly ueidarscood* 

VI. COKCLUDIKC coretnrrs 

This paper haa accacpced co review the key 
tcepa required for cha sodalllng and slaula- 
cion of cha dynaalca and control of future pro- 
poaed large flexible orbiting syste&a which 
will require. In general, both shape aa %«11 
at orlancaclon (acclcudt) concrol. Problaa 
araaa. calnly aasoclaced with cha large order 
of such ayateB codeia, ara highlighted. The 
widespread use of varioua coapucar algorlchas 
required ac dlffarenc acagaa of cha analysla 
should be noted. . 
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Fig. 1 DevelopQcnt of systen software for LSST 
dynamics analysis 
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j-l, p»X, X*-20.52 




— nodAl linoa 
X « frequency par&aator 
j ■ nunber of nodal circles 

p ■ nusber of nodal ceridians 


j-2, p-1, X*-59.86 


Flj* 2 Hade shapes of a free^free circular place (coopucer generated ae Howard Ualverslcy) • 



Fig. 3 Digraph of 16x16 K sa;rlx. 
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III. COOTROLLABILITY OF INHERENTLY DAMPED LARGE FLEXIBLE 
SPACE STRUCTURES 


Abrcract 


I> latTodaetlgti 


' Graph theoretic tccbatquea are txaed to study 
coatrollabUlty of linear systcaa which could re- 
preseat' large flexible orbiting space systeas with 
Inhercat dffjcplng. The coatrollabUlcy of the pair 
of entrlcce representing the systeoi state and con- 
trol Influence catrlccs is assured when states 
In the nsdel are reachable In a digraph sense froa 
at least one Input and also whan Che term rank of 
a Boolean laatrlx whose non trivial cooponenta are 
based on the state and control Influence smtrlces 
baa a tern rank of the order of the state vector. 
It la teen that the damping matrix doea cot In- 
flucsea tha raqulred number of actuate ra but gives 
flexibility to tha poeslbla locations of the actu- 
ators for which tha system la controllabla. and 
Chat Che stiffness oatrlx Cera rank deficiency 
dictates the nusber as veil as Che location of 
the required actuators* Specific example's In- 
clude e model of a shallow spherical orbiting 
shell where both orientation and shape control are 
required, and also a scallac dimensional numerical 
example (unrelated to the shell) which readily 
dcaonstiatea the effect of damping. 


Kopgnelature 


A 

A 







wax aystea matrix 

2 nx2a syetea matrix 

nxB control Influence matrices 

2nxm control Influence matrix 
nxa damping matrlz 
nsn mdlfled damping taatrlx 
Boolean equivalents of A.B end D* 
BBtrlccs 

nxn stiffness matrix 

nxn mass matrix 

reachability matrix 

nxa subsatrlx of matrix R 

diagonal matrix (rxr) 

mxl Input vector 

nxn unitary orthogonal Kitrlx 

nxa unitary orthogonal matrix 

nxl vector 
C3oa matrix 
tcall parameter 

singular value of matrix A 


^Research supported by ?JASA Grant KSC-1414, 
Suppl. 4 

tSenlor Graduate Research Assistant 
*^Prof#ssor of Aefo.'^pace Ejoglneerlng 


Any llAtar. clxse Invariant djnaedcal syscerm 
can be. in general represented by: 

X - AX4BU (1) 

where 

X la en nxl state vector of tha aystea 
A is an nxa systca stsce matrix 
B la an nxa control Influenca matrix 
U la an cod Input of the aystaa. 

sytt» described by equation (1) la aald 
to be contzuXlabla if. with finite U ezxd In finite 
time, the systea (1) can be tranaferrtd froa any 
scat a to any other etate. This concept was first 
Introduced by Richard C. Xalcsn. The verifica- 
tion of coatxoUeblLlty Is essential for control 
systea design as no control law should be designed 
for a syetea which Is not controllable. The eon- 
CzoUablllty concept is even core Important for 
large space structural systems whose dlmensloa- 
allty la very largo. If the design ot the con- 
trol systea Is tadertaken without first verifying 
controllability, a considerable asvuot of effort 
may be wasted, through the fallurr. to arrive et 
any satisfactory control law for an uncontrollable 
eysten. 

In the following sections, the concept of 
controllability, as a property of tha A and B 
matrices, la ravlcwad. 

. .The system. (X) la controllable if and only 

ir*^x 

tank a“"H) - n C2) 

or rank (B.A X^X) » n. !• I.2..*n (3) 

where are the llgcnvalues of the matrix. A« 

The deteminatlon of the rank of the matrix 
In equation (2) paces the problea of selecting n 
Independent colusts out of na colucns and this 
coitld be done by such numerical techniques as the 
singular value decompoaiclon of a taatrlz.^ 

Tha singular value decomposition Is a miaer- 
Icsl algorltha used to find the nuaerlcal rank of 
e reccangular CMtCrlje. A. (axa, say n>a) through 
the evaluation of two orthogonal (unitary) me - 


trices. and such 

that 




(4) 

where p . . -i 

T 


'•[« s] 

• V A7 
1* 2 

(3) 

and S • dlag ' 

.... Op) 
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with b«lng alnsvilar valuta of 

Tbo algea- 

valuaa orthS catrlx AA'^ (or A?A) and r la tha 
rank of tho satrlx A. Tha tachaiffua of finding 
tha alngaUr valuca of A by cval::i\cla8 tha algca- 
valttca of AA^ or A^A ia« la genaral. lota accurate 
than ualag alagular vtlua dacoe>acitlott tach- 
alquas and oay reault in arroncoua concluslona, 
aa la dassnstratod by the following axaapla.^ 



Than [if 1+ti^l* but l^ii^l, for tha cooputar 
accuracy Involved] 


K la the atlffnaaa asatrlx 

Bp la tho control Influence mtrix 

U la tha Qsl vactor of Input a. 


written aa a tat of first order d Iff arantlal aqua- 
tions (la standard state space fora aa): 


-K'^S 


-K-S 


X* r ® " 

X h"^b. 

• . p 


Equation (9) can be considered aa 

X • ax+Fu 


(9) 


( 10 ) 


A^^A ■ 



Where ^ 

(7) * • 


algcnvaluaa of this approxlsatlon to A^A are 
/T and 0, raaptctively, vhareaa the precise 
eigenvalues of a'^a are fo«xd to be and 

U, raepectlvely. 

Tha use of tho singular vmlus decoepoaitloa 
technique, Itsalf.oay rosult In nuaarlcal pro- 
bleas when the aagnltuda of the largtst singular 
value of Che catrl: A, la an order of sa^ltuda 
(or core) different fron thic of the scsallast 
singular value. [A coeplaca discussion of the 
singular value decoeposltloa technique is glvo 
la Ref* 3* It ahould be noted chat with the use 
of this algorlchs, it Is unnecessary to directly 
evaluate the singular values of the A cetrlx. 

The rank deteralnatloa Is eccocpllshcd based on 
the unitary, orthoeoaal properties of the ea- 
trlces Vj_ end Vj within the algorltha.] 



«nd A - D’ - B - H“^B . 

P 

The controllablllcy condlclon (2) for this systea 
can be written es: 

rank fB.Aff,..., - 2n (12) 

If w« aicuaa D-0, which 1» eru* for may Idealized 
frea^^ratlog atructuraa, tha eaoezollabllltr 

C - (b.Ib,..., I Bj, (13) 

bccoaea 


The eppllcaclon of Che conctoUablllcy con- 
dition In equation (3) requires the detentlnatlon 
of Che eigenvalues of the natrlz A and the eval- 
uation of tha rank of a laacrlx of order [ax(a4n)] 
for each of the n eigenvalues* This scheae Is 
oore attractive than thee of condltloa (2) as 
the problea of rank evaluation of an nscia catrlz 
Is reduced to the rank evaluations of n cat rices 
each of diaenslon [nx(a4a)]. The eigenvalues 
ere. In gaierel, ae4Hied for the structural dyna- 
Blc analysis of the systea and cay, thus, be 
already available for this phase of the control 
systoa design* 


II* Controllability of targe Space Structures 

The dynaalcal equations of a large space 
structure systea ere, la geasrel, de:#cribcd by 
a set of linear second order coupled differen- 
tial equatloQ^ as: 

M X DX ♦ rx - Bpff (8) 

where 

X Is Che nxl vector of the generallxed 
coordinates 

M Is the COSO cacrlx (oxn) 

D Is the deep lag matrix (am) (can Include 
viscous daoplug as wall as gyroscopic 
effects) 



It can be very easily se« that C has a rank 2n 
If and only If 

“Ok [B.AB - a (IS) 

vbteh lead* Co ch* foUowla, thMrca^: 

troUable If and only If the pair [A,B] Is con- 
tollable. 

This theorea reduces the decemiaatlon of the 
con cro liability of a 2ath order systea to the de- 
temlnation of the controllability of an equiva- 
lent uCh order systea. la general for large specs 
structure applications, n Itself cay still be suf- 
flccntly large and, thus, nuaarlcal trchnlquos 
would be required in order to determine control- 
lability. This theorea la based on the Inherent 
eseuaption that I>*0 and no Insight can be drawn 
when 0 Is not equal to tero. The effect of the 
Bitrlx. D, on controUability is studied In this 
paper using th-. graph theoretic definition of 
controllability. 
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lit. Graph Theoretic Daflnition 
of Concro liability^ 

Glvea Che g<a«ral, linear, tlioa Invariant 
dy a a nl ca l ayitea described by equation (1), re- 
peated here aa: 

X • AX + BU (16) 

the pair [A.B] la controllable If and otay if: 

(1) the tern rank of 

where are the boolean equlvalenta of the 

satcleaa A and B, respectively; end 

(2) all atatea la the systea are reachable froa 
at Icait one Input in the digraph sense. 

The two tecBSytem rank and reachability, are 
explained here. 



The digraph for the aystea esc rices A, F can be 
drawn in general as ahsvn in Fig. 1. 

the elkaeata of X(i,i4n), -l.e. 

those eleaenta of A appearing In the identity 
ottrix - are represented la the digraph by the 
(solid) lines joining the node^ (i+n) to 1 where 
1-1, 2,..., n. The elements of A(l+a,j), 1-1, 2,.. o, 
j"l»2,...,n l.e. those eleaents of X appearing 
in A ara represented in the digraph by the (dashed) 
lines joining the jth node to tha (I4n)ch node, 
pe el^ta of A(14a,j4n), 1-1,2,.. .,n, j-l,2,..a 
l.e. those elexonts of A appearing In D* are repre- 
sented la the digraph by the (dotted) lines con- 
necting the (j+n)th node to the (14o)th node. 

_ The eleaents of B fthe lover half of F ■ 

B (l^.j)l. Bp(l.j), 1-1,2,. ..n, j-l,2,...a, are 
represented in the digraph by the (double eolld) 
lines joining ths jCh actuator [cho (2n+J)th 
node, hare] to the (l-te)th node. 


Terra rank is the pa xisita rank a eatrix can achieve 
due CO the locations of the non zero, non fixed 
eleaents of the oatrlx rather than dus to the 
nuaerlcal values of the eleaents. A coeplete dis- 
cussion of this concept is provided la the Appei-' 


Beecba blllcy . If one draws a digraph for the ex- 
tended square nstrlx g-i ^xd finds the Input- 

Lo qj 

state reachabUity matrix as explained la the 
Appendix, there cust be at least one non zero 
entry for every row of the submatrlx (R*) la the 
reachability nacrlx, R, forazd from the row and 
column Indices 1 to n, and n+1 to n+o, respectively, 
es shown In equation (17) where n 1m the nuuber 
of states In ths systea and a Is the number of 
actuators. 



(17) 


IV. ControUabUltv of SygretM 
with lahcrtnt Dacrolng 


The dynaalna of large space eysteas with in- 
herent daoplng can be wrlttca as (repeating equa- 
tion (9) with the notatloa defined la equation 


r* 


0 

I 



fo] 

bJ 

m 

A 

D* 


/r 

1*J 


A B 


Froa the reachabUity condition. It ia ob- 
earved that the lines In the digraph due to the 0 
matrix can supplescnt chose lines due to the Bp 
®*^^1** For exasple suppose ciat due to the struc- 
ture of D, there Is a directed (dotted) line fron 
node (n+3) to node (n+2) In Fig. 1. Suppose that 
the actuator represented by node (2n4m) can 
directly laUuence code (tH-3). Then It Is clear 
the actuators number 2 and 3 (represented by node 
2n+2 end 2n+3, respectively) need not be present 
la order to Influence nodes o+2 and n+3. Thus 
Che damping matrix can allow c. greater flexi- 
bility la this selection of the actuator locations. 
Although baaed on this argument cnc could construe 
that the D matrix Influences tfce nuaber of actua- 
tors (that could be removed). It should be reaea- 
bered chat the minima number of ectuators required 
la dictated fron the tern rank, condition which will 
cow be discussed. 

For the system rspresr^ted by eqxiatlon (18) 
to be controllable, the tern rank (aa explained 
la tha Appendix) of the Boolean matrix, m 


0 

I 

— 

0 

s 

“’b 

»E 

0 

0 

0 

- 




mnst be 2n. Kote chat the dlaasislonallcy of the 
state vector ia equation (18) 1m 2a. If A has 
tern rank loss chan n. then tha term rank of the 
(2a-tn)ch order Boolean matrix in equation (19) can 
only be ^gmeaced due to the presence of B», since 
Che Dec lAgJ - - Dec [Ag]. Thtia, Dg can not be 
used to augaant the term rank of the Boolean equi- 
valent of ths state aatrlx A la equation (IS). 

la suissary, the damping marrlx, D, has aa 
effect on the location of the actuators, while the 
matrix A has aa inpact on the Location as well as 
the nunber of actuators. 
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V. yuaerlc^ Eragplei 


The usa cf graph theoretic techaiquaa In the 
detersiinACloa of controllability and the acoent 
of infonsatioa about tha location of tha actuators 
and th« nusibar of actuators noedad Is dosonstratad 
using tha radol of on orbiting shallow spherical 
shall in ccblt with and without tha scabUlxlng 
dumbbell (Fig Z).6 


The Boolean Equivalent of the aystea matrix 
00* for a nhoUov spherical shell is given la 
Fig. 3. The digraph is given la Fig. 4. Froa the 
digraph it can be seen chat c! s nodes may be 
subgrouped as: 


10 

u 

12 

13 

14 

15 16 17 18 

1 

2 

3 

4 

5 

fi 7 8 9 1 


To reach oil the 18 ecatee froa at least one 
Input, control actuators caet directly Influence 
the following nodes: (a) (10 or 11) ; (b) at 

least one of tho nodes. (12-15); (c) (16); (d) (17); 
end (e) (18). The systsn matrix. A. has & tern 
rank deficiency of 1 (note the presence of only 
zeros In the first colusn. Fig* 3) and. thus, one 
actuator Is required for controllability. This 
actuator crust be placed such that the above nen- 
tloned states sre directly Influenced. 


The model of a shallow spherical shell with 
a stabilizing duubbell**' is considered as another 
example for controllability cons tdarac Ions. The 
Boolean €qui\*alcnc of the 22^^ ordered system 
matrix A is given la Fig. 5 and the digiraph is 
shown In Fig. 6. From the digraph it gjtn be seen 
that the total states can be subdivided Into 
two groups 


( 1 ) 


-*•,« ( 2 ). 


12 16 17 18 19 20 21 22 13 14 13 

156789 10 U234 


The control actuators must directly Influence one 
or sore states froa group (1) and one or core 
states froa group (2). Tbs systea matrix here 
has full tera rank and.thu^v one eccuator is suffi- 
cient to establish controllability. 

The two practicO. exssples considered In this 
section to this point do not specifically illus- 
trate tha Independence between the number of 
actuators required and the damping aatrlx. To 
llluetrate this effect an example of sixth order 
is created and analyzed. 


It is asouaed chat Cha ayatea oeeriz la 
given by 


0 0 0 

o 

o 

0 0 0 

0 10 

0 0 0 

0 0 1 

10 0 

1.0 2.0 3.0 

0 0 0 

-10.0 5.0 6.0 

0 0 0 

7.0 8.0 9.0 

. 4 

D J 


The digraph Is drava as shown in Fig. 7 


( 20 ) 
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Froa the digraph the reachability condiclon for 
controllability is satisfied if any one or more 
of the states: 4.5.6. arc directly Infliraccd by 

the control actuators. The term rank of has 
a deficiency of two end thus two actuators are 
required for conctoUabllity. Even if the damping 
matrix 0*0, the same number (2) of actuators is 
needed for controllability and. thus, it is shown 
that damping has no effect on the raqulred number 
of actuators. But. if the damping catrix D«0. 
then the location of the actuators must be cb^ged 
such that the states 4.566 can be dirsecly In- 
fluenced by the control actuators. In order to 
emphasize tha point, when D-0 the dotted lines 
should be resoved froa the digraph slx>wn In Fig* 7. 


VI. Conclusions 

The definition of controllability as 
applied to general linear tine Invariant dynamics 
systems and large space systems is reviewed. 

The special nacura of the coupled matrix second 
order differential equations that are used to 
describe largo space systems is used to airlva st 
specific controllability conditions* The graph 
theory approach is employed to define controllabi- 
lity In terns of the term tank and Input-state 
reachability concepts. This epproach Is used to 
find tha effect of inherent damping present In 
large space systems on the nuaber of tha actuators 
and thoir locations. It is observed that tha 
datsplng does not effect the n-<nfrMa nxmber of 
actuators required, but does provide greater flex- 
ibility In tho possible locations of the actuators. 
Tha number of actuators reqlred depends on the 
term rank of tho generalized system (stiffness) 
matrix. The stiffness matrix also Influences tha 
location of the actuators* 
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Appendix 

A. Tana Rank of a Katrlz 

The cam rank of a square oacrlx of dlsan- 
elon nzn Is laea than o if and only if the eatrix 
has a saro submtrix ”0" of dlscnslon rxr' vlth 
rzr*>n. Tht cam rank is different froa the nu« 
aorlcal rank In tha following seise* If a square 
oatrix of order n has two coluons or rove that ara 
dependdC on each other* then its tara rank ia not 
reduced while its nucsrlcal rank la reduced by one 
for each pair of colums of rows that ara depend* 
enc* For large space systems* the dctarmlnaclDn 
of the aumtrical valuaa of tha alansnts for tha 
system Batrlces ara not exact* and thus the prob- 
ability that tvo colums or rows would ba Iden- 
tically equal* or that ona cow ia a constant tlnaa 
another Is vary stall* If such a dependency exists 
that oust ba detected before subjecting it to the 
term rank tests for establishing tha concrollabi- 
lity of large space systems. 



1 2 3 


Fig. 1. Digraph of [A,B] matrix pair. 


B. Input-State Reachability Matrix 

Tha augmented adjacency catrix for tha ayscaa 
mcrlx pair [A*B] can ba vritten as 


s “ - 


where A^ and ara tha adjacency matrices of 
A and B* respecClTely. 

The states can ba reached from any of tha In- 
puts and can not ba of length core than a* So C« 
can ba raised to the power n and tbua tha augaentad 
syscca reachability isacrix is glvei by 

r 4 1 




where 1* Inpuc-stace reachability matrix 
For all; Che atatas to ba reached from at least one 
Input* every row of R_ oust have at least one non 
rero entry. 
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Fig 5< Location of non-zero eXeocncs of the system matrix of the 
shell with the dumbbell in orbit. 
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IV. EFFECT OF SOLAR RADIATION DISTURBANCE ON A FLEXIBLE BEAM IN ORBIT 
IV. 1 INTRODUCTION 

Proposed future applicatioos of large space structures require control 
of the shape and orientation of the structure in orbit. The principal 
environmental disturbance acting on these structures at the proposed opera- 
tional altitudes are due to the solar radiation pressure. Therefore, it 
is necessary to evaluate the solar radiation pressure effects on the large 
space structures in orbit In order to provide control of their shape and 
orientation. As a specific example of a basic structure, a long flexible 
beam constrained to move only in the orbital plane is considered in this 
study. 

The equations of motion for a long flexible beam oriented along the 

local vertical were obtained previously. Later, the work of Ref. 1 was 

extended to consider the motion and stability of the beam about a nominal 

local horizontal orientation. This system Includes a rigid dumbbell used 

for gravitational stabilization that is connected to the center of mass of 

2 

the beam through a glmballed passive damping device. The control aspects 
of such a beam using point actuators were also considered in Ref. 3. The 
effect of solar radiation pressure on the dynamics of these two types of 
beam structures is studied here, and to the authors* knowledge represents 
the first time that solar disturbance torques acting on large flexible 
Space systems have been treated. 

4 

The force and Ewment expressions obtained by Karymov are used to 
develop the solar radiation disturbance model for a beam by considering 
the individual mode shapes of the free— free beam. The transverse elastic 
displacements are assumed to be small so that the shadowing of the beam 
due to any deflected part of the beam can be neglected. 





OF 


IV. 2 DETERMINATION OF SOLAR RADIATION FORCES AND MOMENTS ACTING ON 
A FLEXIBLE BEAt: 

Let the direction of the incident solar radiation, t, in the body 


coordinate system be denoted as 




(4.1) 


and, let n be the outward unit vector normal to the surface, ds, of a 
body of arbitrary shape exposed to solar radiation (Fig. 4.1). Then, the 
solar radiation force acting on a completely absorbing surface, F^, and that 
acting on a completely reflecting surface,* F^, can be obtained as^. 


F - ht/x-nds 
a o 

s 

F ■-2h / n(t • n) ds 

r o 

8 


(4.2) 


(4.3) 


where, h - 4. 64xl0"* N/m^ is a constant for earth orbiting spacecraft and 
• o 

the integration over an area, s, is bounded by the condition 


T • n > 0 


(4. A) 


The corresponding moments for a completely absorbing surface, N^, and 

— 4 

for a completely reflecting surface, N^, respectively, can be developed as , 


N “ h X X / R(x • n)ds 
a o 

s 

N - 2h / nxR(x • n) ds 
® 8 


(4.5) 


(4.6) 


where R is the position vector of ds with respect to the center of mass. 
For a a surface with an arbitrary reflection coefficient, c^, the force 
and moment expressions become^: 


F - F + e (F -F ) 
er a r r a' 


(4.7) 


N » N + C^(N -N,) 
cr a r r a 


(4.8) 


The forces and moments due to solar radiation pressure acting on a £ree-> 
free flexible beam can now be obtained by considering the shape function 
of the beam, (Fig. 4.1, only the first antl-symmetrlc mode Is depicted). 
The beam is assumed to vibrate In the transverse direction only so that 
the normal at any point is given by 

n - («'l-k)/ /14^ (4.9) 

where ^ and C Is the nondimenslonallzed longitudinal coordinate 

of the beam with the elemental length, 


ds - d£vT+^ ‘ (4.10) 

If the analysis is restricted to a single plane containing ^ and t 
T reduces to 


T ■ a 1+c k 
o o 


(4.11) 


using Eqs. (4.9), (4.10), and (4.11) In Eq. (4.2), the total force acting 
per unit width of the beam is expressed as 


F - -hgT A (a„l+c^k) • (^.'l-k)d5 

a Q O O 

* a^c^h^l+h^c^^k (for symmetric modes) (4.12) 

“ -h^(2a^6 -c )(a 1+c k) for asymmetric modes) 
oooooo 


where, 6 ■ <^!1(0) - deflection at one end of the beam for the n^^ mode. The 
o z 

total force per unit width of the beam acting on a completely reflecting sur- 
face Is obtained after substituting Eqs. (4.9), (4.10) and (4.11) In Eq. (4.3) 
as 




1 


(ao*’-c^)2 

(!+<>• O 


(^’i-k)dC 


(4.13) 
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The expressions for the nooeats per unit width of the beam are developed 
using Eqs. (4.9), (4.10) and (4.11) in Eqs. (4.5) and (4.6) as: 

N - h 7x[ /^(a 4'-c ){(5“7 )i+4ih} dC 

a o ^ o o A 

- -h a c [6 -2/^(fr(Od5] (for symmetric modes) (4.14) 

o o o o Q 

« 2h a c /^«KC)dC (for asymmetric modes) 
o o o o 



2h 

o 

‘•2h 

o 


/V (^0»*-Cq)^ 
o d+(fr'^) 

/I 

o (1-H>^) 


(<fr’i-k)x{(5- |)i+«k) dC 

U'H(5- |)} I dC 


(4.15) 


Eqs. (4.13) and (4.15) Involve complicated line integrals. These Integrals 
can be evaluated using numerical integration methods. For the purpose of 
this numerical study a beam of length 100 meters with tip deflections of 
(1) O.Oli and (ii) O-li were considered. Fig. 4.2 shows the variation of 
the resultant horizontal and normal force components of a beam with a 
completely absorbing surface as the solar incidence angle, 0^^, is varied 
from 0 to 90 degrees. Here, 6^^ represents to angle between the normal 
to the undeflected beam and T. The horizontal and normal force components 
are measured relative to the beam’s undeflected axes. As expected, for 
small tip deflections of the beam, the resultant horizontal absorbing 
force component becomes zero for incidence angles of 0 and 90 degrees, 
respectively, while the normal component has a maximum amplitude at zero 


incidence angle. In Fig. 4.2 and subsequent figures the individual effect 
of each mode, with the assumed beam tip deflection as indicated in the figure, 
is illustrated. 
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J 

f 


'^•'1 
I 

i 

i 

I ^ 

: 4*3 shows the force distribution along the length of the beam deflected 

; In the first nude due to the solar radiation incident at: an angle of 45®. 

\ The asymmetric nature of the force distribution gives rise to a resultant 

% moment about the center of mass of the beam. The magnitude of the resultant 
moments as the solar incidence angle is varied is shown in Fig. 4.4 for each 

\ symmetric mode and the assumed tip deflection. Large cements can result for 

c 

I Xargcr deflections whereas these insnicnts would be zero for a rigid beaiSa 

Because the force distribution for an asymmetric node Is symmetric about 

{ 

j an axis passing through the mass center and parallel tu the incident solar 

radiation, the moments for all asymmetric modes are zero (Fig. 4.3). For 
small pitch angle displacements, the moment due to solar radiation pressure 
may become greater than the moment due to the gravity-gradient forces as 

4 

shown In Fig. A. 5. It is seen that at geosynchronous altitudes, the moment 
due to solar radiation may become predominant even for deflections of the 
order of 0.012.. 

Fign. A. 6, A. 7, and A. 8 show the forces and moments for a completely 
reflecting surface, obtained using numerical integration techniques based on 
Eqs. (A. 3) and (A. 6). It Is seen that the moment for the completely reflecting 
case Increases with the larger value of the tip deflection. Since the radia- 
tion force acts along the normal to the surface for a completely reflecting 
surface, and the deflections of the beam are assumed small, the normal force 
components are seen to be much greater than the horizontal force components 
(Fig. A. 6). Further, the resultant force components also depend on the mode 
shapes (Fig. A. 7) in contrast to the case of the completely absorbing surace 
(Fig. A. 2). 

A. 5 


i 


Bence, the moments for the reflecting beam also depend on the specific mode 
number of the bean Incorporated Into the model as shown in Fig. 4.8. Because 
of symmetric force distribution about the center of mass the resultant moment 
is aero for all asymmetric modes as before. Hence, the moments are zero for 
all asymmetric modes regardless of the surface reflectivity. For the higher 
symmetric modes of the reflecting beam the resultant moments are seen to 
decrease because of the greater scattering associated with sharper changes 
in the beam slope. 

With the aid of these moment diagrams, it is now possible to model the 
disturbance torque due to solar radiation pressure, once the number of modes 
and the associated nodal deflections are specified. This aspect is considered 


In the next section. 


IV. 3 SOLAR RADIATION DISTURBANCE MODEL 


A beam nominally oriented along the local horizontal or local vertical 
is considered. Such a beam makes one revolution per orbit with respect to 
the incident solar radiation. For any symmetric mode and for a given coef- 
ficient of reflectivity, the pitch torque can be expressed as a function 
of the solar incidence angle, 0^^, In the form [from Figs. A. 4 and A. 8], 


where. 


N ■ N sln6. .COS0. 
in 1 X 

N - N + e,(N^ -N^ ) 
D ^ ®m 


(A. 16) 


N ,N “ maximum moment per unit deflection for a completely 

’^m absorbing surface (from Fig. A. 2) and for a completely 
reflecting surface (from Fig. A. 8), respectively. 

For small deflections, N is proportional to the deflection at one end of the 


An(b) 


beam, 6(t), or the nondimens ionalized parameter, e^^ft) - ^ 


, where A^(t) 
’ TL 


modal amplitude function. 

N(t) * e (t) N I sln6, cos0. (4.17) 

n m XX 

and 6^ Is given by 

0.(t) - Cl) t+0(t)+0.(O) (A. 18) 

i' c 1 

where, o)^ Is the orbital angular velocity, and 0 Is the pitch angle of 
the beam. 

The effect of the disturbance on the generic mode is obtained by evaluating 
the integral 

Ejj - e ds (A. 19) 

where, e is the external force. 

Eqs. (4.2) and (A. 3) are substituted Into Eq. (A. 19) to obtain 


E • I • (h^T(fn)} ds 

na z o 

- h c2 

o o o z 


(A. 20) 


A. 7 



and 


- / k*^“^*{-2h^(T*a)2} ds 

- 2h c2 d5 

o o o ^ 


(4.21) 


After combination of Eqs. (4.19) and (4.20), the generic force is obtained 
as. 


E - E +e (E -E ) 

“ °a “r V 

- h c2(1+g) d5 


(4.22) 


For 6j. - 0.5 and a tip deflection of 0.011, Eq. (4.16) yields 

- 2.23xl0"^+0.5(9.4xl0*’^-2.23xl0~^^ 

- 1. 58x10"^ N-m 


This is the maximum torque that is experienced by the beam for a unit de- 


flection equal to lm« in a 100m. length beam at any instant in the orbit. 

The corresponding generic forces on each node with a tip deflection of O.Oli 


In the respective modes are obtained (from Eq. (4.22), for the first four 
niod s, as 

- 0.159x10"® N 

E 2 “ 0.827x10"^ N 

E_j - 0.102x10"® N 
E^ - 0.3432x10"® N 

Thus, the generic forces are seen to be very small and, hence, the nodal 
excitations due to solar radiation pressure are also small. However, the 
magnitude of the solar radiation torques Indicate that considerable pitch 
rotations can be expected for larger deflections of the beam. The numerical 
values for E^ and are used in the two examples in the following sections. 
In the first example a beam nominally oriented along the local vertical 
(Fig. 4.9) is considered. Next, a beam nominally oriented along the local 
horizontal and gravitationally stabilized by a rigid dumbbell (Fig. 4.10) 
is considered. 
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IV.4 EFFECT OF SOLAR RADIATION PRESSURE ON A FLEXIBLE BEAM NOMINALLY 
ORILTTTED ALONG THE LOCAL VERTICAL 


The equations of motion for a thin uniform beam In orbit with its 
axis nominally along the local vertical (Fig* 4.9) is developed In Ref. 2. 

The beam is assumed to undergo only Inplane angular motions and deformations 
and it is assumed also that the center of mass of the beam follows a cir- 
cular orbit. The beam’s elastic motions are considered to be tmconstralned 
and the longitudinal vibrations of the beam are assumed to be negligible 
In comparison with the transverse vibrations. For the case of small amplitude 

pitch oscillations of the beam, the linearized equations of notion are de- 
2 

rived as 


e "+30 




n n n 


En 


M 

n c 


(4.23) 


where. 


0 * pitch motion of the beam 

A 

g «t m non dimen sionalized modal amplitude 
n I 

n * mode number 

J “ pitch moment of Inertia of the beam 

N “ external torque 

E • nth modal force 
n 

M * nth modal mass 
n 

w - orbit angular velocity 
m nodal frequency 



T * oj t, nondimen sionalized time parameter 
c 

t « time 
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Only the first two flexural codes of the beam will be Included in the 
analysis. Using Eqs. (A. 17) and (4.22) and the numerical values for 
and £ 2 , the following three equations of second order result. 

0"+3$ “ 3.6 sin0^ cos0j^ 
e^'+u^e^ “ 3.001x10“^ cos^0j^ 

®2‘*^2®2 “ 1* 563x1 0"** cos^O^ (4.24) 

The second modal oscillation Is seen to be decoupled from the first mode 
and pitch motions. Further, the forcing terms in the first and second modes 
are very small and can be neglected to first order. Therefore, and £2 
have Elutions of the form 

“ c^sinoj^T+c^'costj^T 

^2 * c^sinuy^^+c^costii^^ (A. 25) 

where c^, C 2 » and c^ are constants to be determined from the Initial 

conditions* The pitch equation now becomes 

0**+36 =* 1.8 (c^slncj^T+C 2 ^^s^]['^) sln29^ 

Assuming ®1^^^ ** ^ Q(t) very small 

0.(t) = u t * T from Eq. (4.18) 

X c 

With Cj^(O) - and » 0, Cj^ - 0 and ^2 • c^, and the pitch equation becomes 
0'*+30 » 1.8 p^cjsoj^t sin2T 

“ 0.9e^{sin (2+o)j^)T+sln(2-ai^)T} (4.28) 

The solution of this equation can be obtained in the form 

0(x) - c^sln/3x+CgCos/3T+ slnpx - sinqx 

where, p = 2+ujj^ and q * 2-Uj^ 


4.10 


(4.27) 


With 6(0) ■ 0 and e^(C) « 0.1 and • 10, the pitch response is 
given by 

6(t) - 0. 002392sln/3r-K). 000638slnl2T-0. 001475sln8T 
The response of the beam to the solar radiation disturbance obtained using 
numerical Integration of Eq. (4.24) Is shown In Fig. 4.11* The pitch motion 
shown In Fig. 4.11 Is Identical with the response obtained using Eq. (4.27) 
and shows a maximum pitch amplitude of 0.23®. The effect of the disturbance 
on the first modal oscillations Is seen to be negligible. 
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IV. 5 EFFECT OF SOLAR RADIATIOH PRESSURE ON A DUMBBELL STABILIZED FLEXIBLE 
BEAM NOMINALLY ORIENTED ALONG THE LOCAL HORIZONTAL 

The tacontrolled local horizontal orientation of a bean represents 
an unstable sxstlon. This unstable configuration of the beaa can be 
stabilized by using a rigid duebbell such that the resulting gravity- 
gradient torques provide stabillzatloa. In Ref. 2, the equations of ootion 
for a beam with a dumbbell assumed to be attached at the center of mass of 
the beam (Fig. 4.10) through a spring loaded hinge and having viscous rota- 
tional damping have been developed. In addition to the assumptions made in 
developing Eqs. (4.23), it is further assumed that the duebbell mass is con- 


centrated at the tips and that the viscous force at the hinge is linear. With 
Che usual assumptions of small pitch amplitude and dumbbell oscillations and 
flexural deformations, the linearized equations of motion in Che absence of 


active control and external forces are obtained as^, 

e"->oe’+(k-3)e-ca'-ka+r(7e’-Hcc )C^°^ - , 

n n n z Jai^ 

c 

a"4c,ro'+(c,k+3)a-c,c0’-c,ke-E (ce'H^e )c,C^“^ 

1 X X XnnnXz 

c"+(u2-3)c -S(a-0)4c(o'-9’)}C^“\j /M 2^)+E(c£'-t^G )C^™^ 
nnu zyn — — - 


xa 


m m z 


E /M 
Q n c 


(4.28) 

(4.29) 

(4.30) 


where « J (n,n r l,2y *) and M * mass of the 

z yzz n’’ ’’ n 

beam for all n* 

k * k/J 0 )^; c ■ c/J cx) 
y c y c 

k “ torsional restoring spring constant at the hinge 
c ** viscous damping coefficient 

a « angle between the dumbbell axis and the local vertical 
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I x-0 

beam shape function of the n^^ transverse mode 
Jy/I^ ■ pitch ooment of inertia of the duisbbell 


As before, only the first two laodes will be considered. The forcing terms 

are the sane as for the case of the beam along the local vertical, Eq. (4.24). 
Since the dumbbell is assumed to be rigid, there is no net moment acting on 

the dumbbell due too solar radiation pressure. The first mode influences 
the pitch motion through the forcing function and the second mode affects 
the dumbbell motion through coupling. Thus, pitch, dumbbell, and the two 
modes of the beam are all coupled to each other and the resulting system of 
equations are too complicated to yield analytical solutions. These equations 
were numerically integrated with Initial tip deflections of O.Olt in the 
first mode and zero initial displacements in 0,a and € 2 ^ respectively 
(Fig. 4.12). The steady state response shows pitch amplitudes as high as 2^. 
The first modal oscillations are not greatly affected due to the solar radia- 
tion pressure. The second mode is excited because of the dumbbell motion, 
but the amplitude remains small (maximum * 0.002). The high frequency 

oscillation in ^ pitch acceleration, are suppressed in 

c 

Fig. 4.12 for the sake of simplicity. 

Fig. 4.13 shows the system response for r ctiffer beam with « 20.0 

and the same initial conditions and bean parameters as for the case with 
■ 10.0. The maximum pitch amplitude is seen to be about 0.23^, one 
order of magnitude less than that for the bean with ol « 10.0. 
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In this case the higher frequencies in the second node damp the pitch oscilla- 
tion, through the duisbbell Botlon, inore rapidly so that the pitch amplitudes 
do not build up; Once again, Cj^ and €2 notions are not affected to first 
order because of the solar radiation pressure. Thus, the effect of solar 
radiation pressure Is seen to affect mainly the pitch notion. 

Since, the solar radiation Incidence angle can change considerably 
for synchronous orbits, a long tine simulation (for about 30 orbits) was 
carried out accounting for the change In the Incident angle due to the Earth's 
motion arotuid the sun (-l**/day) as shown In Fig. 4.14. It can be seen that 
errors In both phase and amplitude can result by not Including the annual 
variation in the solar incidence within simulations over long time Intervals. 

The effect of solar radiation pressure on the pitch response for a dif- 
ferent set of Inital conditions (0(0)- a(0)- 0, £j^(0) -62(0)- 0.005) was also 
obtained (Fig. 4.15). The solid line shows the pitch response without the 
®®^®r radiation disturbance. The pitch response in this case Is due to the 
coupled motion In e^, a and 0. Since large amplitude (12°) in pitch motion 
results, the original non— linear equations of motion were used for this study. ^ 
The pitch response in the presence of solar radiation pressure (dashed lines 
In Fig. 4.15) shows a maximum pitch acrolitude of about 9°. Thus, as much 
as 3° difference can result by not including the solar radiation disturbance 
effect for the assumed parameters of the bean In this study. 

Fig» 4.16 shows the effect of solar radiation pressure on a beam which 
is at a low altitude earth orbit (250 n. miles). The pitch excitation is 
seen to be very snail (0.005 ), as expected, because at the low alltudes the 
gravity-gradient torques' are predominant (Fig. 4.5). 
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Fig. 4.6 Variation of Solar Force Conponents with Incidence 
Angle - Totally Reflecting Surface. 
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Fig. 4.9. A Flexible Beau Noninally Oriented Along the Local Vertical 
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Fig. 4.14 Effect of the Annual Variations In Solar Incidence 
Angle- Pitch Response. 
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Fig. 4.16 Pitch Response of Dumbbell Stabilized Flexible 

Bean in the Presence of Solar Radiation Pressure 
at a Low Altitude Orbit. 
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V. HOOP/COLUMN CONTROLS ANALYSES 

V.l DYNAMIC MODEL OF THE HOOP/ COLUMN STRUCTURE 

The structural oodel of the Hoop/Colunn system shown in Fig* 5.1 is 
considered for the control analysis. The hoop is assumed to be constructed 
of five rings with each ring having 24 nodes spaced at 15°. The rings 
are represented by node numbers 1101-1124, 1201-1224, 1301-1324, 1401-1424 
and 1501-1524, respectively. Fig. 5.2 shows a detailed nodal representation 
of the structure including the mast, (nodes 102-127), feeds, (128-136), 
and the solar panels, (99-101). The finite element data of the structure 
provided by the Harris Corporation is used for the controls analysis. 


The dynamic model of the structure can be represented as: 


MX + KX “ B U 
c 


(5.1) 


where X is the state vector containing the generalized coordinates of each 
node and will be of order (nx6) for n number of nodes and all 6 degrees of 
freedom. M is the modal mass matrix of order (6nx6n) and K is the stiffness 
matrix of order '.6nx6n). The control matrix, B„, is of the order of (6nxP) 
for P number of actuators to be arranged on the structure. The data supplied 
by Harris Corporation has eigenvectors for 112 nodes and, therefore, n“112 
for the present model. To decrease the dimensionality of the problem a 
nodal transformation is carried out, by defining 


X - <pq 


(5.2) 


where, is the matrix containing the eigenvectors of Eq. (5.1) and is 
of the order (6nxm) for m number of modes and q is a vector of order (mxl) . 
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Through diagonal Izat ion of Eq. (5.1) the following caatrix eqiiation io 
obtained. 


or 


where 


and 


q + [«’'k 4>1 q - 

r«iJ q + PK^J q - <>\U 


[0^1 


“l 


“2 


m 

m 


[Kj 



(5.3) 


Equation 5.3) is rewritten in the state vector fora as 


'o l‘ 


V 


“ 0 


‘o' 

-K. 



+ 




-i 0 

*1 


• 

q 




u 





— M 


» .. 


The values of K^, and are available with the finite element model. 
The evaluation of the control matrix, for selected actuator locations 
is discussed in the next section. 
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V.2 ARRANGEMENT OF ACTUATORS FOR THE HOOP/COLUKN SYSTEM 

The controls analysis of the Hoop/Colum antenna system requires 
specification of the type of actuators ard their locations and orienta- 
tions In the structure. For this study point thrusters and/or torquers 
are assumed to generate the control forces and torques. The location 
and orientation of these thrusters depend on the mode shapes of the 
structure. The first thirteen nrades corresponding to data provided by 
the Harris Corp. will be Included In the controls analysis and, hence. 

It Is convenient to choose thirteen actuators in the preliminary analysis. 

# 

Each actuator is selected to affect a particular raode, but the same 
actuator Tnay help to control a different node as well* The first six 
nodes are conbinations of rigid body rotations and translations. Actua- 
tors nuicber 5 and 6 are assumed to be arranged as shown in Fig. 5.3 to 
provide control over translation along the x and y directions, respectively, 
and, in addition, also to control the first bending nodes (nodes 8 and >)• 
Actuator 11 controls translation along the z direction, whereas actuators 
8, 12, and 13 control yaw, pitch and roll motions, respectively. Actuators 
1, 2, 3, and 4 are selected so that each actxxator could provide independent 
control of the feed mast torsion (mode 12) • Actuators 9 and 10 are se- 
lected to control the second mast bending (nodes 11 and 13) « Actuator 7 
controls surface torsion (node 10) and is the only actuator assumed to 
be counted on the hoop. The arrangement of these actuators nay need 
reconsideration for more efficient control performance. 



with the selection of in inodes In the model , the dimensions of the 

stSite matrix* A* becomes 2ax2m and the corresponding control matrix* B* 

will be 2mxp for p number of actuators. In the present model, matrix 
I “1 T T T 

B - [0 I [dj^] (> B^] , where <)► (n,d,m) rcspresents the sec of m eigen- 

vectors of Che model and, 

n ■ number of nodes 
d,« degrees of freedom 
m ” number of modes 

^ has a dimension of (112, 6, 13) for the present model and for 13 modes. 
Therefore, the control Influence matrix, Bg, will have a dimension of 
(112, 6, 13) for a total of 13 actuators. The matrix, B^, results from a 
®ls' 26 nt formulation of the load (force and moment) matrix and Is 
developed as follows. A column of the matrix, B^, represents the effect 
of an actuator on the node at which the actuator is located. For example, 
actuator 1 located at node 128 (Fig. 5.3) is assumed to provide a force 
In the y direction only. Pence, the element B^(128, 2, 1) is set equal 
to one and the rest of the elements in the coulumn B^.(nj:, dj, 1) are set 
equal to zero for each n^ and d ^ . The torquer number 8 at node, 98, pro- 
vides only a yaw moment at node, 98, and so the coulumn Bj.(nj^, d^, 8) con- 
tains all zero elements, except at Bj.(98, 6, 8) which is set equal to 1. 
Similarly, the other 11 columns of the influence matrix, B^., are obtained 
as shown in Fig. 5.4 in which the matrix, Bg, is arranged as a two-dimen- 
sional matrix of order (672, 13). Since there are 13 actuators in this 

model, only 13 of the (112x6x13) elements of matrix B are seen to be non 

c 

zero. 

For 13 actuators located as shown in Fig. 5.3 and Table 5.1, the 
T 

given in Table 5.2. The calculation of the was faci- 

litated by the use of a tape containing the elements provided by NASA-LIC. 
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V.3 CONTROLLABILIIY 

control the flntee-eleaent hoop-coluns nodel with 13 iwdee 

1 . th, »del. the alntnom re,p.lred „„a,er of ectnntore is tonm doing 
graph theory. 

Equatioa (5.3) can be cast Into scandard state form as: 




I- 


5.4 


^c^P^r CA.B] in e,„otlon (5.d) Is eontrollohl. If ,„d only If the pnlt 

3 . Cd 1 . eontrollohl,. 1 Fron the reochobllltF condition 
Odd the digraph jhown In Fig. 5.5 for controllability, all the otntea 
oust be Influenced by the Inputs directly. 

The aatrir C ] has a d.ftleeney of 6 In Ito tem rank, aa 

t haa a tern rank of 7. To augnant th, tern rank, (*Tg^] „„„ 

- icaat air linearly Independent non-taro coW, Indleatlng a nlnlno. 
Of am properly plated .etdatora are needed. A poaalbl, aat of actnatora 

are (1,2,3,4.5,12) selected from Table 5.1 or Fig. 53 (v. th 

0 On the contrary, 

Ch. air actnatora ( 1 . 7 . 3 . 4 . 5 . 7 , fron Table 5.1 are not enongh to control 

the thirteen tr>d,a la the ayate. aa at.taa 14,15 and 26 In the digraph 

of fig. 5.5 can not be reached fro. any of the above air l„p„ta („nder 
the aaannptlon that any ,T^^ 

is treated numerically equal to zero) . 


{ ■ 
i ; 
i 


f 

f 
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Actuator no. 

I 2 3 and 4 

5 

6 

7 

8 (torquer) 

9 

10 

II 
12 
13 


Mode being affected 


Feed Mast Torsion (12) 

First Bending (about y axis) (8) 
First Bending (about x axis) (9) 
Surface Torsion (10) 

Taw (rotation about z axis) 
and First Torsion (7) 

Translation along x ^ Also second 
Translation along y Mast bending 
Translation along z 
Fitch (rotation about y axis) 

Roll (rotation about x axis) 


Fig. 5.3 Proposed Arrangement of Actuators - Hoop/Column 
Antenna System 
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j- I VI. DEVELOPMENT OF ALGORITHM TO EVALUATE HOOP/COLUMN COUPLING COEFFICIENTS 

[■' 

[■ 

The generic mode eqiiations and the equations of rotational motion of 
a flexible orbiting body contain coupling terms between the rigid and flexi- 
[ ble modes and terms due to the coupling within the flexible modes that are 

I assumed to be small and, thus, are usually neglected when a finite element 

^ analysis of the dynamics of .the system is undertaken. In this Chapter 

i' a computational algorithm that permits the evaluation of the coefficients 

r in these coupling terms in the equations of motion as applied to a finite 

I element model of the Hoop/Column system is developed. 

\ Using a Newton-Euler approach, one-can express the equations of motion 

I of an elemental mass of the system, in the frame moving with the body, as^ 

I 

^ •• • • 

{a + r+2cii+rHtiKrtiix(w^ =» {f+e4L(q)/p}pdv (6.1) 

f cm 

where p * mass per unit volume, 

e *• external forces per unit mass, 

q “ elastic transverse displacements of the 
element of volume. 

f “ force due to the gravity on the unit mass, and 

L « the linear operator which when applied to q yields 
the elastic forces acting on the element of volume 
considered. 

r * position vector of element dv 

(0 ■ inertial angular velocity of the body frame 
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VI. 1 EQUATIONS OF ROTATIOlUL MOTION 


The equations of rotational isotlon of the body are obtained by taking 


the moments of all the external. Internal and Inertial forces acting on 
the body, l.e., from Eq. (6.1) 


/ r X [a + r + 2ancr + (oxr + ux(ajxr)]pdv 
cm 


/ rx[L(q)/p+ f + ejpdv 


( 6 . 2 ) 


one can obtain the following form for the equations of rotational 


motion. 


R + Z + Z ^ 

n“l n=l n«l 


(6.3) 


where R “ /[r x(ujxr ) - (r * 0 )) (user )] pdv 
Y o o o o 

•• # • • • 

f Q > /[r xq+2r x(o)xq)+T x(o)xq)+qx(cjxr^) 
n=l V ® 


-(r^*o)) (a3xq)“(q*ti)) (o)xr^)] pdv 
L » /q pdvx(a -f ) + ? to^A / r x $ 

n-i ^ n-1 ^ ^ V ° 


G„ * / r xMr pdv 
R o o 


? G “ /[r xMq+qxMr ] pdv 


C * /rxe pdv 


r - r^+q ^ 

M - matrix operator'*’ which when applied to r yields gravity-gradient forces 

a _ « acceleration of the center of mass 
cm 

* force/mass due to gravity at the undeformed center of mass 
^ (n)« modal shape vector for the n^^ mode 

0)^ * frequency of the n^^ mode 

A * time dependent modal amplitude function 
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VI. 2 GENERIC MODE EQUATIONS 


The generic code e<iuacioQ Is obcslned by csking the DodeX components 
of ^1 1 internal, external and Inertial forces acting on the body, i.e. , 


/ ♦ • [a +r+2uncrhjjxr+uDc(ujxr) ] pdv 

V 

» [L (q)/p)+f+el pdv 


(6.4) 


The generic mode equation is obtained in the following form 5 

A +u*A +4^/H + I “ Cg^+ - 

n n n n n n a nn n n . n 

where Cf - ‘uxr + ❖ *iI»c((lKr )] pdv 

n Y o o 

£ (C > •UKq+^ *t«q+^ 'toxCukq)] pdv 

m-1 no V 

R, - ‘Mr^pdv; Z g “ /T^”^*Mqpdv; 

Ti o uj*}^ xan Y 

E « /T^”^*epdv and D' - /'^^"^pdvCa -f ). 
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VI. 3 CARTESIAN COMPONENTS OF . TOE DIEFEBENT COUPLING TERMS 

^ S r ^ 

The expressions for R» Q r ^ n, 

in Cartesian components are presented in this 


^ tan' ^n' ^mn 
section. 


One can express the following vectors in 
their Cartesian component form as 

.. r^“C:ti+5y3+5sk? « 


n«l 


iind 


5(n) 


g(n) 


Q^n>i ^ + Q 


(n): 


.(n): 


G'"'i * * 01 "’h. 


where i, j , k are unit vectors along the body principal axes 
of inertia in the undeforaed state; Cy» Cg 

ordinates of a point in the undeformed state. 

With the use of the component forms of the 
vectors given above, one can expand the various vector e-v- 
pressions given in Bqs. (6.3) and (6.5) to obtain 


+ tJyWy + <’^3c'^2^"z“x^ ^ 
(J.Uj + t Jy-^x’ 


( 6 . 6 ) 
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' zx xz '”z ‘“y'"z* zz yy ' “’x"’y' xz . 

+ u (u 

zx ' X z jcy yx ■ '^z "y* ' yz 


+H^"' ) 1 


(6.7) 




( 6 . 8 ) 




+ M3^(H<“>+H^>) + 2Mj3 W^' -H^;' ) I 


(6.9) 


a X yz zy y zx xz ^ z xy yx 

+ M a + a a ) •+ a a 

X y' xy yx ' y z yz zy ' z x zx 

xz ' X yy zz ^ y' zz XX ' 


"z' XX yy ’ 


(6.10) 


<D - 2A (a + a (L<°"> 

^mn x' yz zy * y' zx xz ' 


+ u (l"“'-l"“')] + a [i (l'^’-l*™”) 

z xy yx m X 'yz zy ' 

• + i +a 

y' zx xz ' z xy yx ’ 
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where and 

(n) 


(n) 


way when a is y in H^g I*, 


a, B » X, y/ z or 1/ 2 , 3. t^hen a is x in H^g «*. «^g 

the corresponding value of o in M^g is 1* In a si m ilar 

(mn) 


or L 


(mn) 


'oB 


o is z in H^g* or « is 3 in The same rea- 

soning holds for 3 also. 


, a is 2 in M^g and when 

A 

aB* 


The expressions for and are 

obtained by the cyclic permutation of x, y» z in the 
expression for in Eq- (6.7) and the expressions for 


,{n) 

^y 


(n) 

^z 


and G2“' are obtained by the cyclic permutation of 
X/ yy z in the expression for in Bq. (6.9). 
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For a discretized model the expressions for the volume integrals are 
replaced by the following summations* 

( 6 . 12 ) 

' tJ '1 1 

(a, 6 - x,y,z) 
k 




■-a'r’' 


1-1 “ ^ 


(n), 

6 'i“l 


where 


k - total number of discrete masses 
1 - index identifying a nodal point 

m^^ - mass concentrated at the i^^ node. 

^ - coordinates of m^^ in the undeformed state 


( 6 . 13 ) 
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VI. 4 EVALUATION OF COUPLING COEFFICIENTS IN THE EQUATIONS OF MTION AS 
APPLIED TO A FINITE ELEMENT WDEL OF THE HOOP/COLUMN SYSTEM 

VI. 4.1 Model Description 

The structural dynamic modeling of the Hoop/Column antenna has gone 
through many stages before reaching the single surface model which will 
be analyzed in this chapter. 

Initially, it had 231 nodes distributed as follows: 192 nodes on 

the 8 support circles including the hoop (24 nodes on each circle spaced 
at 15^ intervals); 28 nodes on the mast and the feed mast; and 11 nodes 
at the points of location of the solar panels (upper and lower), the 
S band reflector, and the feed panels (up-llnk and down-link)»*8ee Figs. 
5.1 and 5.2. After reduction the number of nodes was diminished to 114 
including a total of 96 nodes on the circles: 1100, 1200, 1300, and 1400 
7 nodes on the mast and the feed mast; and 11 nodes at the locations of 
the solar panels, the S band reflector, and the feed panels (Fig. 6.1). 


VI. 4.2 Approximate Mass Distribution 

2 

From an unpublished document prepared by the Harris Corporatlont 
and submitted by NASA Langley Leuearch Center, it has been possible to 

I 

arrive at the mass distribution shown in Table 6*1. 9803*0 lb. out of | 

the total weight of the Hoop/Column Antenna (10»070 lb.) were distributed 

b<!tween the final grid points. The distribution was done in agreement with 

the information found in the Harris Corporation document. The page numbers 

appearing in Table 6.1 refer to particular mass/moment of inertia calcula- 

. 2 

tions in the Harris Corporation document. 

The small (2%) discrepancy between the calculated total mass (9803.0 lb.) 
and the stated weight of the system (10,070 lb) is thought to be attributed 
to: (1) uncertainties in the weight of specific stringers; (2) uncertainties 

inherent with the finite element reduction technique where the initial mass 
must be redistributed between a reduced, final number of grid (node) points; 
and (3) other miscellaneous uncertainties, such as the exact weight/location \ 

I 

of the optical instrument, etc. \ 

VI. 4.3 Cartesian Coordinates of all the Nodal Points in the Final 
NASTRAN Output 

Reference 2 contains the cylindrical coordinates of all the nodal points 
on the mast, the feed mast, and at the location of the panels and electronics. 

It also contains the Z coordinates of the planes which contain the circles 
along with their respective diameters. Thus, the Cartesian coordinates of 

all the nodal points were obtained by a simple transformation from cylindrical i 

! 

to Cartesian coordinates. i 
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VI. 4. 4 Development of a Coiqputational Algorithm for Evaluation of 
the Coupling Coefficients 

After receipt of the tape containing the modal functions, this 
information was stored in our IBM 360 in such a manner that when one 
calls subroutine^GETMP(2) , he can refer to the component of the 
mode shape vector at the grid point J by VECMP(I,J,K). Based on this, 
an algorithm described in the flow diagram, Fig. 6.2, was designed and 
tested. As indicated in Fig. 6.2, the available data, such as: the 

Cartesian coordinates of the grid points on the mast, the feed mast and 
Che ones at the locations on the appendages; and such as the mass concen- 
trations at all the nodal points are input into the software routine and 
these data will consequently have to be updated according to any develop- 
ment in the Hoop/Column modeling. The subroutine, DCS, (given the radius 
of the circles and the Z component of their centers) computes the Cartesian 
coordinates of the nodal points on the circles. 

Subroutine GETMP(2), which makes the available^ is called and 

the values of components of the desired mode shape vector at the particular 
grid point are Incorporated into a loop mathematically described by Eqs. 
(6.12) and (6.13). It should be noted that, for reasons of effectiveness, 
each coefficient is evaluated separately on the circles and on the other 
grid points and then combined to yield the corresponding coupling coefficient 
for the entire Hoop/Column system. 

The algorithm has been tested for two modes (the 7^^ and the 8^^) 
successfully, but only after the evaluation of the coefficients corresponding 
to all the 13 modes will one be able to make positive conclusions. 
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Given X,Y,Z (cartesian coordinates of the nodal 
points on the mast and feed) and the mass concen- 
tration at each nodal point in the system (approxi- 
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Hoop assembly at grid points 1101, 1107, 1113, TOTAL 9803 O 

1119; 244.51bs/point ' 

Table 6.1 Approximate Mass Distribution at Final Grid Points (Pounds) 











































VII. GENERAL CONCLUSIONS AND RECOMMENDATIONS 

The widespread use of various computer algorithms required at 
different stages for the simulation of the dynamics and control of 
large flexible orbiting systems should be emphasized. Problem areas 
are mainly associated with the large order required to model such systems. 
The use of graph theoretic techniques can often be used to reduce the 
computational effort involved in the calculation of the eigenvalues of 
such large ordered systems. Computer generated interactive graphics 
can provide additional insight into the interpretation of the flexible 
modal shape functions of complex systems. 

The graph theory approach can also be utilized to define controll- 
ability in terms of the term rank and input-state reachability concepts. 
This approach can be employed to examine the effects of inherent damping 
(usually expected to be present in LSST systems) on the number and loca- 
tions of the required actuators. It is seen that the damping matrix does 
not Influence the required number of actuators but offers greater flexi- 
bility to the possible locations of the actuators for which the system is 
controllable. The system (stiffness) matrix term rank deficiency dictates 
the number of actuators required and also influences the location of the 
actuators. 

A mathematical model of the solar radiation forces and moments 
acting on a free-free flexible beam in orbit has been developed. For 
small pitch angles, it is seen that the solar radiation torques due to 
the deformations of the beam can be larger than those due to the grcvity— 
gradient for orbits near synchronous altitude. 
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In^orbix-plane steady state open- loop (uncontrolled) responses for different 
initial beam deflections Indicate that, in general, the effects of solar 
pressure on the modal amplitudes are small, but the magnitude of the 
induced pitch oscillations can be relatively larger. Future work could 
extend the model to plate and shell surfaces, and also assess the effect 
of the solar pressure disturbance on previously developed control lows, 
designed primarily to provide certain transient response characteristics. 

A preliminary analysis of the finite element dynamic model using 
the first 13 inodes of the 122m. Hoop/Column antenna system indicates 
that a minimum of six properly placed actuators is required for controll- 
ability. Additional work is currently underway to analyze transient 
responses and force-impulse requirements for control laws based on dif- 
ferent techniques using ORACLS, and also various cc binations of number 
and location of actuators. 

Finally, an algorithm has been developed to evaluate the various 
coupling terms between the rigid rotational and flexible modes and also 
the intra-modal coupling terms in the equations of motion using the Hoop/ 
Column mass distribution as provided by the NASTRAN finite element pro- 
gram as an example. Such coupling terms are usually not included in 
finite element models which are based on the earth-based vibratic.ial and 
rigid modes only. Current work is in progress to evaluate the order of 
magnitudes of these terms - at least for the first 13 modes that are being 
used in Hoop/Column controls analysis. 
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